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Abstract 


Composite  materials  are  currently  being  considered  for 
nvimerous  Naval  applications.  Among  their  many 
advantages,  their  inherent  vibration  damping 
characteristics  is  becoming  a  dominant  design 
consideration.  To  date,  however,  there  are  no  analytical 
techniques  which  have  the  capability  of  predicting  the 
damping  loss  factor  of  a  particular  composite  system. 
One  is  therefore  limited  to  experimental  determination 
of  the  damping  loss  factor  for  composite  materials  which 
requires  a  high  degree  of  accuracy  and  precision.  An 
experimental  apparatus  has  been  designed,  fabricated  and 
tested  using  2024  T-4  aluminum  as  a  calibration 
material.  This  material  has  been  well  characterized  in 
regard  to  its  damping  loss  factor.  Of  critical 
importance  is  the  material  placed  between  the  clamping 
block  and  the  specimen.  These  materials  act  to 
distribute  the  load  and  to  minimize  frictional  losses. 
Various  materials  were  investigated.  The  material  that 
provided  results  which  were  the  most  consistent  and 
matched  those  analytically  anticipated  was  a  fiberglass 
woven  supported  teflon  sheet.  Experimental  results  for 
the  2024  T-4  aluminum  calibration  material  match  those 
determined  analytically  to  within  5%  over  a  frequency 
range  of  40  to  300  Hz.  The  variation  in  results  at  a 
particular  frequency  are  within  10%. 
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INTRODUCTION 

There  have  been  numerous  investigators  who  have  developed 
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test  techniques  for  the  determination  of  the  damping  loss  factor 
for  composites.  (1-5)  For  the  experimental  determination  of  the 
damping  loss  factor  of  composites,  the  most  recently  utilized 
techniques  are  the  forced  resonant  vibration  technique (1,2)  and 
the  impulse  technique. (3 , 4 , 5)  Of  these  techniques,  the  impulse 
technique  appears  to  be  the  most  widely  used.  In  the  impulse 
technique,  the  specimen  is  normally  held  in  a  cantilever  beam 
configuration,  and  is  excited  via  an  impact  near  the  clamped  end. 
The  beam  response  is  measured  using  either  accelerometers  or  a 
noncontact  eddy  current  probe.  From  the  beam  displacement  versus 
time  information,  the  frequency  response  function  is  determined. 
Using  the  half  power  band  width  method,  the  loss  factor  is  then 
calculated. 

The  manner  in  which  the  sample  is  supported  is  of  vital 
importance  to  the  accurate  determination  of  the  material  damping 
loss  factor.  Fixturing  and  materials  that  are  utilized  should 
dissipate  a  minimal  amount  of  energy.  Loss  that  can  occur  at  the 
clamping  area  include  frictional  losses  and  dissipation  by  the 
material  used  to  clamp  the  specimen. 

There  are  numerous  technical  reasons  why  the  impulse 
technique  has  gained  wide  acceptance  for  testing  composites. 
First,  all  vibrational  frequencies  are  excited  when  the  specimen 
is  subjected  to  an  impact.  Numerous  fundamental  frequencies  and 
the  various  modes  of  vibration  can  therefore  be  investigated 
using  a  single  specimen.  This  is  accomplished  by  varying  the  beam 
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length.  Secondly,  to  minimize  external  energy  losses  from  such 
sources  as  friction  at  the  clamped  region  or  from  aerodynamic 
damping,  beam  tip  displacements  must  be  minimized.  For  low 
frequencies,  which  are  of  major  concern  for  various  structural 
applications,  placing  the  beam  into  resonant  forced  vibration 
results  in  large  tip  displacements.  The  damping  loss  factor  that 
is  therefore  determined  is  a  combination  of  the  material  internal 
damping  and  external  losses.  This  information  would  therefore  be 
of  little  use  for  design  purposes.  With  an  impact  excitation, 
minimal  force  levels  can  still  be  used  with  resultant  tip 
displacements  being  kept  to  a  minimum. 

This  paper  describes  the  experimental  apparatus  developed  for 
the  experimental  determination  of  the  damping  loss  factor  for 
composite  materials,  including  both  the  fixturing  and  the 
hardware  and  software  used  for  data  acquisition.  In  addition,  the 
results  of  the  experimental  calibration  of  this  system  is 
reported,  where  2024  T-4  aluminum  beams  are  used  as  the 
calibration  specimen.  This  material  system  was  chosen  since  it 
has  been  experimentally  well  characterized  and  has  been  shown  to 
analytically  follow  Zener  thermoelastic  theory (6)  for  the  damping 
loss  factor  as  a  function  of  frequency.  (7,8) 

The  experimentally  determined  damping  loss  factor  results  are 
then  compared  with  the  analytically  determined  loss  factor  using 
Zener  thermoelastic  theory. 
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EXPERIMENTAL  APPARATUS 


For  vibration  damping  testing,  there  are  two  primary 
considerations  when  designing  fixturing  for  testing  materials. 
First,  it  is  necessary  that  the  specimen  be  isolated  from  its 
surroundings.  No  vibrational  energy  from  external  sources  should 
be  allowed  to  influence  the  vibrational  response  of  the  specimen 
being  tested.  Accomplishment  of  this  likewise  infers  that  the 
vibrational  energy  imparted  to  the  specimen  will  not  be 
dissipated  by  the  fixturing  due  to  an  energy  transfer  from 
the  specimen.  Secondly,  care  must  be  taken  to  minimize  all 
other  possible  sources  of  energy  dissipation  so  that  the 
measured  damping  is  the  material  inherent  damping  loss 
factor.  Two  prevalent  sources  for  external  energy  dissipation 
which  should  be  minimized  are  frictional  losses  at  the 
clamping  regions(2)  and  aerodynamic  damping. (7 , 9-11) 

In  order  to  isolate  the  specimen  from  the  surroundings,  the 
specimen  fixturing  was  attached  to  the  ] oad  frame  of  a  Sontag 
fatigue  test  machine,  model  SF-IU.  This  machine  consists  of  a 
36  X  36  X  3  in.  solid  steel  top  attachment  plate  to  which  is 
attached  the  components  for  fatigue  testing.  The  entire 
system  is  isolated  from  the  machine  frame  via  four  springs. 

The  entire  system,  therefore  is  isolated  from  the 
surroundings.  In  addition,  the  mass  of  this  support  system  is 
several  hundred  pounds,  which  will  be  not  be  easily  excited 
by  the  vibration  of  the  test  specimens. 
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The  fixturing  that  is  used  for  the  testing  is  shown  in  figure 
1.  This  consists  of  a  steel  base  plate,  12  x  15  x  1  in.,  which  is 
bolted  to  the  top  plate  of  the  Sontag  machine  with  four  0.75  in. 
diameter  bolts.  A  solid  aluminum  breadboard  is  then  bolted  to 
this  plate  with  four  screws.  The  breadboard  consists  of  a  12  x  12 
inch  plate  with  1/4-20  mounting  holes  on  a  1  inch  spaced  sguare 
pattern  starting  1.5  inches  from  the  edge.  A4x3xl.25in. 
clamping  block  base  is  then  attached  to  the  breadboard  with 
three  screws.  This  clamping  block  base  was  machined  to  accept  two 
0.75  in.  diameter  tool  steel  guide  posts.  These  are  used  to 
ensure  that  the  specimen  will  be  clamped  flat.  On  top  of  the 
clamping  block  base  is  positioned  the  clamping  block.  This 
consists  of  a  4  x  1.25  x  1.25  in.  steel  plate  which  is  machined 
to  accept  linear  bearings  through  which  the  guide  posts  will 
pass.  Test  specimens  with  widths  up  to  1.5  in.  and  of  any 
thickness  can  be  clamped  by  this  fixture  using  two  1/4-20  screws. 

The  eddy  current  probe  is  positioned  in  proximity  to  the  beam 
tip  using  the  X-Z  vernier.  The  vernier  is  attached  to  the 
breadboard  with  four  cap  screws.  This  specific  vernier  has  a  0.5 
inch  travel  in  each  direction.  By  repositioning  the  vernier  at 
different  bolt  locations  on  the  breadboard,  beams  of  any  length 
can  be  interrogated  with  the  eddy  current  probe.  It  should  be 
noted  that  a  Plexiglas  fixture,  to  which  the  eddy  current  probe 
is  bolted,  was  fabricated  to  attach  the  eddy  current  probe  to  the 
vernier.  The  use  of  Plexiglas  was  required  since  the  eddy 
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current  probe  was  unshielded  and  metallic  structures  near  the 
sensor  can  result  in  spurious  displacement  determination.  A 
schematic  of  the  entire  assembly  is  given  in  figure  2. 


EXPERIMENTAL  PROCEDURE 

For  damping  loss  factor  determination,  beam  tip  displacement 
versus  time  information  is  required.  The  eddy  current  probe  is 
used  to  measure  the  tip  displacement.  This  probe  is  electrically 
excited  and  outputs  a  voltage  dependent  on  its  distance  from  a 
conducting  medium.  This  particular  probe  can  be  used  for  both 
magnetic  and  non-magnetic  materials.  The  probe  has  a  measuring 
range  of  0.120  inches  with  a  resolution  of  0.01%  of  the  measuring 
range.  The  probe  must  first  be  calibrated  before  use.  To 
accomplish  this,  it  is  placed  into  a  micrometer  calibration 
fixture,  Kaman  part  number  850854-001.  Following  the  calibration 
procedure,  the  output  of  the  eddy  current  probe  is  adjusted  to 
the  requirements  of  the  test.  The  gain  can  be  adjusted  to 
optimize  sensitivity  for  particular  tests.  For  this  test,  the 
probe  was  calibrated  so  that  1  volt  was  equal  to  0.040  in.  of 
displacement.  Displacement  vs.  voltage  information  was  obtained 
using  the  fixture  in  increments  of  0.005  in.  A  linear  fit  was 
performed  on  the  data  and  the  slope  obtained  for  use  in  the  data 
reduction  programs. 

To  conduct  a  test,  a  specimen  is  first  placed  into  the  test 
fixture.  The  eddy  current  probe  is  positioned  so  that  its  center 
in  approximately  1  in.  from  the  specimen  end  and  in  the  middle  of 
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the  specimen  width. 


The  output  from  the  eddy  current  probe  and  the  force  hammer 
are  attached  to  the  MetraByte  DAS-16  screw  terminal  accessory 
board  using  BNC  connectors.  The  accessory  board  is  connected  to 
the  DAS-16  data  acquisition  board  installed  in  an  IBM  PC-AT.  This 
is  a  16  channel  high  speed  A/D  interface  with  direct  memory 
access  (DMA) .  This  system  has  a  10  MHz  on-board  crystal 
controlled  oscillator  to  provide  the  time  base.  The  system  can  be 
set  up  to  read  +  5  volts.  The  information  is  stored  and  read  in 
bits  with  a  range  of  -2048  to  2047  bits  corresponding  to  -5  and 
+5  volts,  respectively.  The  system  resolution,  then,  with  the 
eddy  current  probe  calibration  set  to  0.04  in.  per  volt,  is 
9.78  X  10  in. 

A  computer  program  was  written  in  basic  which  controls  the 
DAS-16  data  acquisition.  The  program  utilizes  numerous  call 
routines  which  are  provided  with  the  DAS-16  board.  The  program 
initially  lets  the  user  adjust  the  height  of  the  probe  so  that 
the  output  of  the  probe  is  zero.  The  program  then  prompts  the 
user  for  the  appropriate  attachments  and  acquisition  rate.  The 
ability  to  vary  the  acquisition  rate  allows  the  user  to  zoom  into 
a  particular  frequency  range.  In  order  to  prevent  aliasing,  or 
distortion  of  the  Fourier  transform(12) ,  however,  the  acquisition 
rate  should  be  at  least  8  times  the  frequency  of  interest. 

The  data  acquisition  system  takes  the  desired  number  of 
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data  points  at  the  prescribed  rate  which  can  be  varied  up  to  a 
maximum  of  76  kHz.  The  data  that  is  taken  is  stored  in  a  specific 
memory  location  which  can  be  accessed  to  perform  a  screen  plot  of 
the  log  of  the  displacement  versus  time  and/or  displacement 
versus  time.  In  the  former  case,  the  plot  allows  the  user  to  make 
a  quick  check  on  the  amplitude  decay  to  ensure  qualitatively  that 
there  are  no  secondary  sources  of  energy  dissipation.  This  would 
be  indicated  by  the  appearance  of  two  distinct  slopes  of  the  peak 
height  values  from  the  log  displacement  versus  time  plot.  Figure 
4  is  an  example  of  the  log  displacement  versus  time  plot  from  the 
vibration  response  of  2024  T-4  aluminum  subjected  to  an  impact 
excitation.  Figure  5  is  the  displacement  versus  time  plot  of  the 
same  data  sample. 

The  user  can  likewise  continue  through  the  computer  routine 
without  plotting  the  data.  At  this  point,  the  program  prompts  the 
user  for  the  file  name  and  address  where  the  data  will  be 
transferred.  The  data  can  be  transferred  as  the  digital  stored 
information,  or  it  can  be  stored  as  actual  displacements.  An 
option  in  the  program  allows  the  user  to  input  the  specific 
scaling  factors  obtained  from  the  calibration  of  the  eddy  current 
probe.  The  program  then  creates  and  opens  the  storage  file  and 
transfers  the  appropriate  data  there.  In  addition,  the  program 
places  a  header  in  the  data  file  that  includes  the  acquisition 
rate  and  interval  between  data  points.  A  annotated  copy  of  this 
data  acquisition  program  is  given  in  Appendix  1. 
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To  perform  the  Fast  Fourier  Transform  (FFT)  on  the  data,  the 
Interactive  Laboratory  System  from  Signal  Technologies 
Incorporated  is  utilized.  The  system  is  first  used  to  import  the 
data  as  a  sample  file.  The  system  partitions  this  data  into 
records  of  up  to  1024  data  points.  The  system  can  then  be  used  to 
perform  the  FFT  on  any  record  and  the  information  is  stored  to  a 
file.  The  system  is  then  used  to  transfer  the  FFT  information  to 
an  ASCII  file.  The  information  that  is  transferred  includes  the 
frequency,  dB  magnitude  and  phase  for  each  data  point. 

This  file  is  then  read  into  a  word  processing  software  called 
Wordstar.  The  frequency,  dB  magnitude  and  phase  values  near  the 
resonant  frequency  value  are  read  into  two  files.  One  file 
contains  approximately  twenty  data  points  in  which  the 
frequencies  are  less  than  the  resonant  frequency,  and  the  second 
file  contains  the  twenty  data  points  which  are  greater  than  the 
resonant  frequency.  The  values  for  each  of  the  two  data  sets  can 
be  easily  distinguished  from  the  phase  information.  In  FFT 
analysis,  there  is  a  2  tt  phase  shift  at  the  frequency  which 

corresponds  to  a  resonant  frequency.  A  sample  of  this  phase  shift 
is  given  in  figure  6,  where  the  dB  magnitude  versus  frequency  and 
phase  angle  versus  frequency  are  given. 

The  two  data  sets  are  then  used  as  input  to  a  basic  program. 
The  program  converts  the  dB  magnitude  information  to  log  values. 
These  converted  files  are  then  read  into  a  graphics  routine  by 
Golden  Software,  Inc.  called  Grapher.  Each  data  set  is  fit  with  a 
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fourth  order  polynomial  using  orthogonal  coefficients.  The 
orthogonal  coefficients  for  the  data  set  given  in  figure  6  are 
shown  in  table  1. 

The  specific  fourth  order  equation  with  orthogonal 
coefficients  that  is  used  for  the  curve  fit  routine  is  given  as 
follows 

y  =  B, +(x-a2  )  [Bj +(x-a3  )  (B3 +(x-a4  )  (B^ +(x-a5  )B5  )  -  B5]  - 

p3(B,+[X  -tts  ]B5  )]  -  P3  [B3+(X  -a,  ){B,  +(X  -as  )B5  )  -P^  Bj]  (1) 


where  the  B,  are  the  orthogonal  coefficients,  and  the  a,  and  p^ 


are  the  recursion  factors.  The  value  of  x  given  in  equation  1  is 
calculated  from  the  users  x  value  of  the  frequencies  using  the 
following  equation 


=  r .  ^  ^ 

l^[i.e.  fr< 


input 


-  X. 


frequency]  "”‘<1  I 


(2) 


where  the  x^^j^  and  are  scaling  factors  for  the  original 
frequency  data. 


A  file  is  then  created  with  the  values  of  the  orthogonal 
factors,  and  the  alpha  and  beta  recursion  factors.  This  file  is 
then  used  as  input  to  a  basic  prograirf  which  calculated  the 
intersection  of  the  two  curves,  i.e.  the  left  and  right  curves 
for  the  resonant  frequency  peak,  and  the  frequencies 
corresponding  to  the  half  power  points  for  the  specific 
curves.  The  half  power  points  are  the  values  of  the  frequency 
that  correspond  to  .707  of  the  magnitude  of  the  peak  height  value 
for  the  resonant  peak.  The  loss  factor  is  then  determined  using 
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the  following  equation 


A  copy  of  the  computer  program  which  calculates  the  loss 
factor  using  the  fourth  order  polynomial  fit  with  orthogonal 
coefficients  is  given  in  Appendix  2. 


EXPERIMENTAL  RESULTS 

In  order  to  calibrate  the  system  and  to  determine  the 
accuracy  and  precision  of  the  procedure,  2024  T-4  aluminum  was 
tested.  Beams  were  machined  from  0.125  in.  plate  to  specimens 
with  widths  of  1.0  inches  and  lengths  of  12  inches.  The  12  inch 
lengths  were  machined  parallel  to  within  0.001  in.  to  minimize 
any  errors  in  the  test  results. 


Beam  lengths  were  determined  that  would  result  in  various 
first  resonant  frequencies  using  the  following  equation (9) 


N 


Elg 
W,  L‘ 


X  c, 


(4) 


where  E  is  the  bending  stiffness,  I  is  the  moment  of  inertia  of 
the  beam,  g  is  the  acceleration  due  to  gravity,  W,  is  the  weight 

per  unit  length,  L  is  the  beam  length  and  C,  is  the  constant 
for  the  i'^’mede  of  vibration.  (13) 


As  was  previously  mentioned,  there  are  two  sources  of  energy 
dissipation  that  must  be  minimized  to  ensure  that  the  damping 
loss  factor  experimentally  determined  is  of  the  material  internal 
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damping  and  not  from  external  dissipation  sources.  To  minimize 
aerodynamic  damping,  loss  factors  were  only  determined  for 
samples  where  maximum  tip  displacements  were  less  than  +  0.008  in. 
This  represents  a  displacement  less  than  7%  of  the  total  beam 
thickness.  This  was  chosen  since  it  was  previously  shown  that 
aerodynamic  damping  could  be  minimized  if  tip  displacements  were 
kept  less  than  20%  of  the  specimen  thickness. (10, 11) 

The  second  source  of  energy  dissipation  that  can  be 
controlled  and  should  be  minimized  are  the  frictional  losses  at 
the  clamped  region  of  the  beam.  An  investigation  was  conducted 
where  various  materials  were  placed  between  the  specimen  and  the 
clamping  block.  These  materials  included  tapered  aluminum 
clamping  blocks,  tapered  fiberglass  blocks,  fiberglass  fabric 
reinforced  polyurethane  sheet,  fiberglass  fabric  reinforced 
teflon,  and  no  clamping  block.  These  systems  were  positioned 
between  the  sample  to  be  tested  and  the  clamping  blocks.  The 
specimens  were  tested  using  the  procedure  described  above.  The 
experimentally  determined  loss  factors  were  compared  to  that 
analytically  anticipated  from  Zener  Thermoelastic  Theory.  The 
clamping  material  which  resulted  in  the  lowest  experimental  value 
of  the  loss  factor  for  the  beam  was  the  fiberglass  fabric 
reinforced  teflon  sheet.  The  loss  factor  determined  using  this 
material  was  also  the  best  match  to  that  analytically 
anticipated.  Since  all  other  clamping  materials  gave  results  that 
were  higher  than  anticipated,  it  was  assumed  that  they  were 
providing  additional  source  for  energy  dissipation.  The  material 
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that  was  used  for  the  rest  of  the  testing  was  therefore  the  fiber 
glass  fedsric  reinforced  teflon  sheet. 

Beam  lengths  were  varied  so  that  the  first  resonant  frequency 
was  varied  from  approximately  38  Hz.  to  300  Hz.  The  test 
procedure  given  above  was  followed  and  the  loss  factors 
determined.  The  results  from  the  testing  are  given  in  Table  2. 


The  results  of  the  test  program  were  compared  with  the 
predictions  determined  using  Zener  thermoelastic  theory . (1, 2 , 4 , 7) 
The  2024  T-4  aluminum  has  been  shown  to  follow  predictions  for 
the  loss  factor  given  by  this  theory.  The  loss  factor  is  given  by 
the  following 

w  T 


_  ot^  E  T  _ 

^  (1  +  t2  ) 

where  a  is  the  coefficient  of  thermal  expansion,  E  is  the 


(5) 


bending  stiffness  of  the  beam,  T  is  the  absolute  temperature,  C 
is  the  heat  capacity,  w  is  the  angular  frequency  and  t  is  the 

relaxation  time  for  the  heat  flow  across  a  beam  of  rectangular 
cross  section.  The  value  of  t  is  determined  by 


T  = 


h^  c 

7t2  K 


(6) 


where  h  is  the  beam  thickness  and  K  is  the  thermal  conductivity 
of  the  beam  material. 


For  each  of  the  beam  lengths  used  in  the  testing,  the  first 
resonant  frequency  was  analytically  determined  using  equation  4. 
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These  results  are  compared  with  those  determined  experimentally 
in  Table  2.  All  of  the  experimental  values  of  the  first  resonant 
frequency  are  within  the  error  for  the  analytically  determined 
values . 

Also  given  in  Table  2  are  the  experimentally  determined  loss 
factors.  For  comparison  purposes,  the  anticipated  loss  factor  of 
the  specimens  using  Zener  thermoelastic  theory  are  given.  Table  3 
gives  the  values  used  for  input  to  equations  5  and  6  for  the 
determination  of  the  damping  loss  factor  for  the  aluminum 
calibration  specimen.  The  experimental  loss  factors  are  within 
2.5%  of  those  analytically  anticipated  for  frequencies  less  than 
196  Hz.  At  frequencies  greater  than  196  Hz.,  the  experimental 
results  were  within  10%  of  those  anticipated  analytically. 

The  results  from  the  testing  are  plotted  along  with  the 
analytical  curve  of  the  loss  factor  from  Zener  thermoelastic 
theory  in  Figure  7.  In  this  figure,  the  solid  curve  is  the 
analytical  values  of  the  loss  factor.  A  fourth  order  curve  fit 
was  performed  on  the  experimental  data  to  determine  the  degree  of 
fit  of  the  testing  at  frequencies  between  those  at  which  tests 
were  conducted.  A  plot  of  the  curve  fit,  also  given  in  figure  7 
shown  by  the  dotted  line,  shows  excellent  agreement  with  the  plot 
of  the  analytically  determined  values  of  the  loss  factor.  As  can 
be  seen  from  the  plot,  the  curve  fit  for  the  data  matches  the 
results  predicted  by  Zener  thermoelastic  theory  to  within  2%  in 
the  range  for  which  data  was  taken. 


CONCLUSIONS 


The  experimental  apparatus  that  was  designed  for 
determination  of  the  damping  loss  factor  for  materials  in  the 
cantilever  beam  configuration  accurately  determines  the 
damping  loss  factor  with  a  high  degree  of  precision.  Resonant 
frequencies  determined  from  curve  fitting  of  the  frequency 
response  function  are  within  the  error  of  the  analytically 
determined  resonant  frequencies.  The  variation  in  discrete  values 
of  the  loss  factor  at  a  specific  frequency  have  variations  of 
less  than  14%  with  the  exception  of  the  data  at  312  Hz.  where  the 
spread  was  22%.  The  anticipated  trends  in  loss  factor  for  the  2024 
T-4  aluminum  are  determined  to  within  2%  through  appropriate 
curve  fitting  of  the  data  compared  to  Zener  thermoelastic  theory 
for  frequencies  from  38  to  100  Hz,  and  less  than  5%  for 
frequencies  between  100  and  312  Hz.  Secondary  sources  of  energy 
dissipation  have  therefore  been  minimized. 
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Table  1:  Fourth  order  polynomial  fit  statistics  for  peak  values 
from  the  FFT  given  in  figure  6. 


LEFT  SIDE  OF  PEAK 


Degree 


Orthogonal 

Factors 


Recursion  Factors 
Alpha  Beta 


0 

1 

2 

3 

4 


-3.6147 

1.24339 

0.458577 

0.161536 

0.0916427 


0 

-0.000321833 

0.000065269 

-0.000107365 

0.000147984 


0 

1.55585 

1.22238 

1.1425 

1.07872 


RIGHT  SIDE  OF  PEAK 


Degree 


Orthogonal 

Factors 


Recursion  Factors 
Alpha  Beta 


0 

1 

2 

3 

4 


-3.86925 

-1.39188 

0.496194 

-0.234782 

0.122776 


0 

-0.000219959 

-0.000141179 

0.000124755 

0.000279605 


0 

1.52375 

1.20287 

1.13346 

1.08324 
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Table  2:  Experimentally  determined  vibration  damping  loss  factor 
results  for  2024  T-4  aluminvim  tested  in  cantilever  beam 
configuration. 


Calculated 

Resonant 

Frequency 

(Hz) 


Experimental 

Resonant 

Frequency 

(Hz) 


Experimentally 
Determined 
Loss  Factor 
(  X  10  “*  ) 


Average 
Loss  Factor 

(  X  lO'**  ) 


Loss  Factor 
from  Zener 
Thermoelastic 
Theory 
(  X  10  “’  ) 


38.85+1.40 


48.91+1.76 


78.29+2.82 


80.38+2.89 


96.82+3.49 


195.20+7.03 


199.73+7.19 


223.92+8.06 


312.08+11.24 


38.563082 

38.57834 

38.58244 

49.0575 

49.0637 

49.0639 

49.0684 

49.0706 

49.0909 

49.5670 

77.17365 

77.18895 

77.20165 

77.2040 

77.25243 

80.1034 

80.1891 

80.2092 

80.2496 

95.328 

96.47 

197.2326 

197.2516 

197.1784 

198.1793 

198.2263 

198.5398 

219.628 

219.645 

219.856 

220.041 

305.1664 

305.1771 

305.1800 

305.1902 


19.66 

21.91 

25.45 

22.12 

21.99 
22.79 

20.29 

24.85 
23.02 
22.57 

20.30 

21.24 
18.52 

20.85 

20.25 

20.00 

21.28 

19.99 

20.88 

18.44 

18.60 

11.85 
11.35 
10.08 

9.84 

8.62 

10.31 

8.28 

8.59 

9.35 

10.09 

8.00 

7.69 

5.83 

8.80 


22.34 


22.74 


20.23 


20.54 


18.52 


11.09 


9.59 


9.08 


7.58 


22.21 


22.94 


20.70 


20.04 


18.62 


10.90 


10.64 


9.61 


7.05 
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Table  3:  Input  values  used  in  equations  3  and  4  for  the 

determination  of  the  damping  loss  factor  for  the 
aluminum  calibration  specimen  using  Zener  thermoelastic 
theory . 


INPUT 

PARMETER 

COEFFICIENT 
OF  THERMAL 
EXPANSION 

BENDING 

STIFFNESS 

TEMPERATURE 

(ABSOLUTE) 

HEAT 

CAPACITY 

THERMAL 

CONDUCTIVITY 

SYMBOL 

a 

E 

T 

C 

K 

VALUE 

23.2x10 

73.097 

298.72 

2560511.2 

126.47 

fin/ in 

GPa 

(°K) 

r  J  1 

r  kg  m  ^ 

1  “K  J 

Im'  °  K  J 

Lsec  ^  °  K  J 
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Figure  1:  Schematic  diagram  of  experimental  apparatus  for  vibration  damping  testing. 


Figure  2:  Experimental  apparatus  used  for  the  vibration 
damping  testing  with  2024  T-4  Aluminum  specimen. 
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using  lictsie  luograin 


2024  1  -4  Al.  (400  Hz) 


Log  displacement  versus  time  information  for  2024  T- 
Aluminum  sample  data 
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h-MASt  (RADlAliS)  MAGniTUDE  (oB) 


Figure  6:  FFT  of  2024  T-4  Aluminum  showing  dB  magnitude  versus 
frequency  and  phase  angle  versus  frequency. 
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2024  T-4  ALUMINUM  LOSS  FACTOR 


FREQUENCY  (Hz) 

;  factor  for  2024  T-4  Aluminum  determined  using  Zener 
rrnoelastic  Theory  and  from  experimental  testing. 
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APPENDIX  1 


:  An  annotated  Basic  computer  routine  used  for  the 
data  acquistion  system  to  read  displacement 
versus  time. 
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10  '  THIS  PROGRAM  IS  THE  DATA  AQUISITION  PROGRAM  FOR  THE  VIBRATION  DAMPING 
20  '  PROJECT.  THIS  IS  SET  UP  TO  READ  CHANNEL  0  AS  THE  DISPLACEMENT  SENSOR 

30  •  AND  CHANNEL  1  AS  THE  FORCE  HAMMER.  THE  PROGRAM  TAKES  DATA  AT  YOUR 

40  '  PRESCRIBED  RATE  AND  WRITES  IT  TO  A  DATA  ARRAY,  CH%  AND  DT.  THE  PROGRAM 
50  '  PLOTS  THE  ABSOLUTE  VALUE  OF  THE  LOG  OF  THE  DISPLACEMENT  VS.  TIME.  THE 
60  '  DATA  CAN  BE  SAVED  IN  EITHER  THE  DIGITAL  VALUE  FORM  (-2048  TO  2047)  OR 

70  '  CAN  BE  SCALED  BEFORE  BEING  SAVED  TO  A  FILE.  THE  DIGITAL  DATA  FILE  NAME 

80  '  WILL  BE  VIBDATA.DAT  WHILE  THE  SCALED  DATA  WILL  BE  VIBDATA.PRN.  THE  PROGRAM 
90  '  WILL  ALLOW  YOU  TO  MAKE  THE  APPROPRIATE  SCALING  FACTOR  CHANGES. 

100  '  ********************************************************************** 

110  '  THIS  PROGRAM  IS  CALLED  VIBDAMP.BAS 
120  '  ********************************************************************** 

125  DIM  DT% (12000) 

126  DIM  CH% (12000) 

127  KEYPRESSED  =  0 

128  CLS 
130  " 

140  '  LOAD  DASH16.BIN  DRIVER 
150  ' 

160  DEF  SEG  =  &H5000 

170  BLOAD  ''C:\ILS\DAS\DASH16.BIN",  0 
180  • 

190  '  ***  INITIALIZE  WITH  MODE  0 
200  " 

210  MD%  =  0 

220  DIO%(0)  =  &H330  •  SET  I/O  ADDRESS 

230  DI0%(1)  =  2  ^  SET  INTERRUPT  LEVEL 

240  DIO%(2)  =  3  •  SET  DMA  LEVEL 

250  FLAG%  =  0 
260  DASH16  =  0 

270  CALL  DASH16  (MD% ,  D10%(0),  FLAG%) 

280  IF  FLAG%<>0  THEN  PRINT  "INSTALLATION  ERROR"; STOP 
290  ' 

300  '  ***  SET  MULTIPLEXER  SCAN  LIMITS 
310  MD%  =  1 

320  'INPUT  "ENTER  LOWER  SCAN  LIMIT,  (0-7)  : " , LL% 

330  'D10%(0)  =  LL% 

340  'INPUT  "ENTER  UPPER  SCAN  LIMIT,  (0-7)  UL% 

350  'DI0%(1)  =  UL% 

360  PRINT 

370  PRINT  "  ATTACH  THE  OUTPUT  FROM  THE  EDDY  CURRENT  PROBE  TO  CHANNEL  0" 

380  PRINT  "  AND  THE  OUTPUT  FROM  THE  FORCE  HAMMER  TO  CHANNEL  1" 

390  INPUT  "  AFTER  DOING  THIS,  HIT  ANY  KEY  ",Z$ 

395  CLS 

400  D10%(0)  =  0 
410  DI0%(1)  =  0 

420  CALL  DASH16  (MD%,  D10%(0),  FLAG%) 

430  IF  FLAG%<>0  THEN  PRINT  "ERROR  #  " ;FLAG% ; "IN  SETTING  SCAN  LIMITS": STOP 
440  ' 

450  '  ***  DO  ONE  A/D  CONVERSION  AND  INCREMENT  MUX 

460  ' 

470  PRINT 
480  'MD%  =  3 

490  'CALL  DASH16  (MD%,  D10%(0),  FLAG%) 

500  'PRINT  USING  "CHANNEL  ##  READS  #####  BITS" ; DIO% ( 1) ,  D10%(0) 

510  'INPUT  "DO  YOU  WANT  TO  DO  ANOTHER  ?",  A$ 

520  'IF  A$  =  "Y"  OR  A$  =»  "y"  THEN  410 

530  '  ***  DO  1266  A/D  CONVERSIONS  AND  PRINT  AVERAGE 

540  PRINT 

550  PRINT  "  THE  FOLLOWING  SECTIONS  WILL  ALLOW  YOU  TO  ADJUST  ZERO  FOR  THE" 
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560  PRINT  "  EDDY  CURRENT  PROBE.  ADJUST  THE  VERNIER  FOR  THE  PROBE  HOLDER" 

570  PRINT  ’*  UNTIL  THE  OUTPUT  IS  CLOSE  TO  0.  THEN  HIT  ANY  KEY  TO  CONTINUE" 

580  WHILE  NOT  KEYPRESSED 
590  MD%  *  17 

600  D10%(0)  =  10  '  DIVIDE  10  MHz  BY  10  TO  GIVE  IMHz  FREQUENCY 
610  DI0%(1)  =  13  '  DIVIDE  1  MHz  BY  12  TO  GIVE  83.3  KHz  FREQUENCY 
620  FLAG%  =  0 

630  CALL  DASH16  (MD%,  DIO%(0),  FLAG%) 

640  D10%(0)  =  1266 

650  DI0%(1)  =  &H6800  •  SEGMENT  OF  MEMORY  TO  RECIEVE  DATA 

660  DIO%(2)  =  1  'INDICATES  TYPE  OF  TRIGGER;  1  =  TIMER,  0  =  EXTERNAL  TRIGGER 
670  DIO%(3)  =0  '  0  *  ONE  CYCLE,  1  =  RECYCLE 

680  MD%  =  6 

690  CALL  DASH16  (MD%,  DIO%(0),  FLAG%) 

700  MD%  *  8 

710  CALL  DASH16  (MD%,  D10%(0),  FLAG%) 

720  IF  DI0%(1)  *  1  THEN  GOTO  700 
730  MD%  =  9 

740  D10%(0)  =  1266  '  NUMBER  OF  WORDS  TO  TRANSFER 

750  DI0%(1)  =  &H6800  '  SEGMENT  OF  MEMORY  TO  TRANSFER  FROM 

760  DIO%(2)  =  0  '  SET  TRANSFER  TO  BEGIN  AT  BEGINNING  OF  SEGMENT 

770  DIO%(3)  =  VARPTR(DT%(0) )  '  TO  START  TRANSFER  AT  BEGINNING  OF  ARRAY 

780  DIO%(4)  =  VARPTR(CH%(0) )  '  CHANNEL  THAT  DATA  IS  FROM 

785  CALL  DASH16  (MD%,  DIO%(0),  FLAG%) 

790  IF  KEYPRESSED$<>""  THEN  KEYPRESSED=-1 

800  SUM  *  0 

810  FOR  I  =  1  TO  1266 

820  SUM  =  SUM  +DT%(I) 

830  NEXT  I 

840  ZERO®  SUM/1266 

850  LOCATE  15,20 

860  PRINT  USING  "###.###"; ZERO 

870  KEYPRESSED$®"" 

880  KEYPRESSED$=INKEY$ 

890  WEND 

895  ' 

896  '  ***  SET  MULTIPLEXER  SCAN  LIMITS 

897  MD%  =  1 

898  D10%(0)  =  0 

899  DI0%(1)  =  1 

900  CALL  DASH16  (MD%,  D10%(0),  FLAG%) 

901  IF  FLAG%<>0  THEN  PRINT  "ERROR  #  " ;FLAG% ; "IN  SETTING  SCAN  LIMITS"; STOP 

909  ' 

910  '  ***  SET  PROGRAMMABLE  TIMER  RATE  USING  MODE  17 

920  '  SO  THAT  TIMER  CAN  BE  SET  TO  MAX.  VALUE 

930  '  NOTE  THAT  CURRENT  TIME  FOR  2  CHANNEL  SCAN  IS  1.30  E-5  SEC." 

940  ' 

950  PRINT  "DO  YOU  WANT  TO  CHANGE  THE  VALUES  FOR  THE  PROGRAMMABLE  TIMER" 

960  PRINT  "  THE  CURRENT  SETTING  IS  FOR  A  FREQUENCY  OF  76,9  KHz.  THIS  IS" 

970  PRINT  "  THE  MAXIMUM  FREQUENCY.  ANY  NEW  FREQUENCY  WILL  BE  LOWER  THAN" 

980  INPUT  "  THIS.  DO  YOU  WANT  TO  CHANGE  FREQUENCIES  ";A$ 

990  PRINT 

1000  IF  A$  ®  "N"  OR  A$  ®  "n"  THEN  GOTO  1140 

1010  PRINT  "  TO  SET  A  NEW  FREQUENCY,  INPUT  TWO  NUMBERS  WHICH  WILL" 

1020  PRINT  "  BE  DIVIDED  INTO  10  MHz  THAT  WILL  RESULT  IN  THE  FREQUENCY" 

1030  INPUT  "  THAT  YOU  DESIRE;  A,B  (CURRENT  DIVISOR  IS  130)  ";D,V 

1040  K  »  D*V 

1050  PRINT 

1060  TI#-  1E+07/K 

1070  IF  K>-130  GOTO  1170 
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1080 

1090 

1100 

1110 

1120 

1130 

1140 

1150 

1160 

1170 

1180 

1190 

1200 

1210 

1220 

1230 

1240 

1250 

1270 

1280 

1290 

1300 

1310 

1320 

1330 

1340 

1350 

1360 

1370 

1380 

1390 

1400 

1410 

1420 

1430 

1440 

1450 

1460 

1470 

1480 

1490 

1500 

1510 

1520 

1530 

1540 

1550 

1560 

1570 

1580 

1590 

1600 

1610 

1620 

1630 

1640 

1650 

1660 

1670 

1680 


PRINT 

PRINT  "  YOU  MUST  CHOOSE  A  DIVIDE  BY  GREATER  THAN  130  IN  ORDER  TO” 
PRINT  "  COLLECT  DATA  CORRECTLY.  PLEASE  INPUT  VALUES  WHICH  HAVE  A" 
PRINT  "  PRODUCT  GREATER  THAN  130  !  !  !  '• 

PRINT 
GOTO  1010 
D»5 

V  =  26 

TI#=  lE+07/130 
MD%  =  17 

DIO%(0)  =  D  •  DIVIDE  10  MHz  BY  10  TO  GIVE  IMHz  FREQUENCY 
DIO%(l)  =  V  '  DIVIDE  1  MHz  BY  12  TO  GIVE  83.3  KHz  FREQUENCY 
FLAG%  =  0 

CALL  DASH16  (MD%,  D10%(0),  FLAG%) 

/ 

#  ***  DO  N  A/D  CONVERSIONS  AND  TRANSFER  TO  MEMORY  VIA  DMA  -  MODE  6 


INPUT  "ENTER  NUMBER  OF  CONVERSIONS  DESIRED  (UP  TO  12000):  ",M 
PRINT 

PRINT  "THE  SYSTEM  WILL  BEGIN  TAKING  DATA  FROM  ASSIGNED  CHANNELS" 

INPUT  "WHEN  READY,  STRIKE  ANY  KEY  AND  DATA  COLLECTION  WILL  BEGIN  ",A$ 
D10%{0)  =  M 

DI0%(1)  =  &H6800  •  SEGMENT  OF  MEMORY  TO  RECIEVE  DATA 

DIO%(2)  =  1  'INDICATES  TYPE  OF  TRIGGER;  1  =  TIMER,  0  =  EXTERNAL  TRIGGER 
DIO%(3)  =0  '  0  =  ONE  CYCLE,  1  =  RECYCLE 

MD%  =  6 

CALL  DASH16  (MD%,  D10%(0),  FLAG%) 

MD%  =  8 

CALL  DASH16  (MD%,  DIO%(0),  FLAG%) 

IF  DI0%(1)  =  1  THEN  GOTO  1360 


PRINT  CHR$(7) 


'***  TRANSFER  DATA  FROM  MEMORY  TO  ARRAY  USING  MODE  9 
MD%  =  9 
D10%(0) 

DI0%(1) 

DIO%(2) 

DIO%(3) 

DIO%(4) 


CALL  DASH 16  (MD%, 
PRINT  CHR$(2) 


M 

&H6800 

0 

VARPTR(DT%(0) ) 
VARPTR(CH%(0) ) 
D10%(0) 


NUMBER  OF  WORDS  TO  TRANSFER 
SEGMENT  OF  MEMORY  TO  TRANSFER  FROM 
SET  TRANSFER  TO  BEGIN  AT  BEGINNING  OF  SEGMENT 
TO  START  TRANSFER  AT  BEGINNING  OF  ARRAY 
CHANNEL  THAT  DATA  IS  FROM 
FLAG%) 


'  READ  DATA  FROM  MEMORY  SEGMENT  AND  PRINT  GRAPHICS,  LOG (DISPLACEMENT) 
'  VS.  TIME  AND  FORCE  VS.  TIME 


FOR  I  =  1  TO  M-1  STEP  2  '  SEARCH  HAMMER  ARRAY  FOR  WHEN  VALUE 

IF  DT%(I)  <  250  THEN  GOTO  1600  '  REACHES  A  PRESET  VOLTAGE  INDICATING 

B  =  I  -  21  '  HAMMER  IMPACT.  STORE  ARRAY  NUMBER  AS  B 

K  =  (M-B)/640  '  K  SCALES  TIME  TO  FIT  ON  X-AXIS 

GOTO  1620 
NEXT  I 
B  »  M-21 
Q=0 

FOR  I  =  0  TO  M  STEP  2  '  SEARCH  FOR  MAXIMUM  VALUE  OF  DISPLACEMENT 

IF  ABS(DT%(I))  <  Q  THEN  GOTO  1660  '  SO  THAT  GRAPHICS  CAN  BE  SCALED 

Q  »  DT%(I) 

NEXT  I 
R#-Q/50 

INPUT  "DO  YOU  WANT  TO  PLOT  LOG (DISPLACEMENT)  VS.  TIME  ";A$ 
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1690  IF  A$  =  ''N”  OR  A$  =  "n”  THEN  GOTO  1820 
1700  CLS: SCREEN  2: KEY  OFF 

1710  W  =  (LOG(Q) )/100  '  W  SCALES  DISPLACEMENT  TO  FIT  ON  Y-AXIS 

1720  FOR  I  =  B  TO  M-1  STEP  2  •  LINES  1230  -  1410  PLOTS  DISPLACEMENT 

1730  Z  =  ABS(DT%(I))  '  AND  FORCE  VS.  TIME  USING  LOG  PLOT 

1740  IF  Z  *  0  THEN  GOTO  1780 
1750  L  =  LOG(Z) 

1760  PSET  ((I-B)/K,  (100  -  (L/W) ) ) 

1770  GOTO  1790 

1780  PSET  ((I-B)/K,  100) 

1790  PSET  ((I-B)/K,  (190- (DT% (I+l)/23) ) ) 

1800  NEXT  I 

1810  INPUT  '"',A$ 

1820  INPUT  "DO  YOU  WANT  TO  PLOT  THE  SAME  DATA  USING  DISPLACEMENT  VS.  TIME";A$ 

1830  IF  A$  =  "N"  OR  A$  =  "n"  THEN  GOTO  1920 

1840  CLS  .'SCREEN  2  .'KEY  OFF 

1850  FOR  I  =  B  TO  M-1  STEP  2 

1860  PSET  ((I-B)/K,  (51- (DT% (I)/R#) ) ) 

1870  PSET  ((I-B)/K,  51) 

1880  PSET  ((I-B)/K,  ( 195- (DT% ( I+l) /23 ) ) ) 

1890  NEXT  I 

1900  N  =  M-B 
1910  INPUT  •"',A$ 

1920  PRINT  "NUMBER  OF  DATA  POINTS  IS  ";N 

1930  INPUT  "WOULD  YOU  LIKE  TO  CHANGE  THE  START  OF  THE  DATA  ";A$ 

1940  IF  A$  =  "N"  OR  A$  =  "n"  THEN  GOTO  1970 

1950  PRINT  "THE  CURRENT  VALUE  OF  THE  START  FOR  DATA  IS  ",'B 

1960  INPUT  "  ENTER  NEW  VALUE  FOR  START  OF  DATA  TRANSFER  ";B 

1970  CLS; SCREEN  0 

1980  • 

1990  •***  TRANSFER  DATA 
2000  ' 

2010  PRINT  "  THE  DATA  WILL  NOW  BE  PRINTED  TO  A  FILE.  INSERT  A  DATA  DISK  IN" 
2020  PRINT  "  DRIVE  A.  INPUT  THE  NAME  YOU  WISH  TO  CALL  THIS  FILE" 

2030  PRINT  "  PLEASE  USE  THE  EXTENSION  OF  'PRN'  FOR  YOUR  FILE 

2040  INPUT  "",FILE$ 

2042  B$  =  ":" 

2044  IF  INSTR(FILE$,B$)  =  0  THEN  GOTO  2047 

2045  C$=FILE$ 

2046  GOTO  2050 

2047  A$="A:" 

2048  C$=A$+FILE$ 

2050  PRINT 

2060  PRINT  "  CURRENT  ARRAY  WILL  BE  TRANSFERRED  TO  A  FILE,  ";C$ 

2065  PRINT 

2070  PRINT  "  YOU  CAN  SCALE  THE  DATA  OR  PRINT  THE  DIGITAL  DATA  TO  THIS  FILE" 
2075  PRINT 

2080  PRINT  "  THE  SCALING  FACTORS  ARE  CURRENTLY  SET  WITH  A  " 

2090  PRINT  "  SLOPE  AND  INTERCEPT  FOR  DISPLACEMENT  OF  3.989063E-2  AND  " 

2100  PRINT  "2.897717E-5,  RESPECTIVELY.  DO  YOU  WANT  TO  PRINT  THE  DIGITAL" 

2130  INPUT  "DATA  TO  FILE";A$ 

2140  IF  A$="Y"  OR  A$="y"  THEN  GOTO  2340 
2145  PRINT 

2150  PRINT  "DO  YOU  WANT  TO  CHANGE  THE  SLOPE  AND  INTERCEPT  OF  THE  EDDY" 

2160  INPUT  "CURRENT  PROBE  FROM  3.989063E-2  AND  2.897717E-5,  RESPECTIVELY" ;A$ 
2170  IF  A$  *  "N"  OR  A$  =  "n"  THEN  GOTO  2210 
2180  PRINT 

2190  INPUT  "  INPUT  NEW  SLOPE  AND  INTERCEPT  IN  FORMAT,  S,I  ";S#,E# 

2200  GOTO  2230 

2210  S#»  3.989063E-02 
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2220  E#=  2.89771E-05 

2230  OPEN  C$  FOR  OUTPUT  AS  #1 

2240  PRINT  #1/'THE  FREQUENCY  OF  TEST  IS  ••;TI#/2 

2250  PRINT  #1,"THE  TIME  INTERVAL  BETWEEN  DATA  POINTS  IS  ";2/TI# 
2260  FOR  I  =  B  TO  M  STEP  2 
2270  DT#  =  DT%(I)/409.4 
2280  DT#  =  S#  *  DT#  +  E# 

2290  PRINT  #1,DT# 

2300  NEXT  I 
2310  CLOSE  #1 

2320  PRINT  "THE  NUMBER  OF  CONVERSIONS  THAT  HAVE  BEEN  MADE  IS  ";N 
2330  GOTO  2435 

2340  PRINT  "THE  DIGITAL  DATA  WILL  BE  STORED  IN  FILE  ";C$ 

2350  PRINT  "THE  NUMBER  OF  DATA  POINTS  ARE  ";M-B 
2360  OPEN  C$  FOR  OUTPUT  AS  #1 

2370  PRINT  "THE  FREQUENCY  OF  DATA  ACQUSITION  IS  ";TI#/2 
2380  FOR  I  *  B  TO  M-1  STEP  2 
2400  PRINT  #1,  DT%(I) 

2410  'PRINT  #1,T#;DT%(I) ;DT%(I+1) 

2420  NEXT  I 
2430  CLOSE  #1 
2435  PRINT 

2440  PRINT  "WOULD  YOU  LIKE  TO  RUN  THIS  PROGRAM  AGAIN  ? 

2450  INPUT  "",M$ 

2460  IF  M$  =  "Y"  OR  M$  =  "y"  THEN  GOTO  127 
2470  END 
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APPENDIX  2 


:  Basic  computer  routine  used  for  the  determination 
of  the  resonant  frequency  and  damping  loss  factor 
from  the  fourth  order  curve  fit  values  using  the 
half  power  band  width  method. 
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10  '  ******  the  following  program  is  called  40RTH0FIT  ***** 

20  ' 

30  '  THIS  PROGRAM  PERFORMS  A  FIT  TO  DATA  IN  A  DATA  FILE  FOR  A 
40  '  FOURTH  ORDER  POLYNOMIAL  FIT  WITH  ORTHOGONAL  COEFFICIENTS 
50  '  THE  INFORMATION  IS  READ  FROM  A  FILE  AND  YOU  ARE  PROMPTED  FOR 
60  •  THE  APPROPRIATE  INFORMATION 
70  ' 

80  PRINT  "  IN  ORDER  TO  RUN  THIS  PROGRAM,  YOU  NEED  TO  HAVE  A  DATA  FILE  " 

90  PRINT  "  THAT  HAS  THE  FOURTH  ORDER  POLYNOMIAL  COEFFICIENTS  ALONG  WITH" 

100  PRINT  "  THE  ALPHA  AND  BETA  VALUES.  THEY  SHOULD  BE  INPUT  FOR  EACH" 

110  PRINT  "  EQUATION  IN  THE  FORM:  " 

120  PRINT  "  A,  B,  C,  D,  E" 

130  PRINT  "  Al,  A2,  A3,  A4 ,  A5  " 

140  PRINT  "  Bl,  B2,  B3 ,  B4 ,  B5" 

150  PRINT 

160  PRINT  "  WHERE  THE  FORM  OF  THE  EQUATION  IS  EX^4+DX^3+CX^2+BX+A  " 

170  PRINT  "  AND  Al  -  A5  AND  Bl  -  B5  ARE  THE  ALPHA  AND  BETA  COEFFICIENTS" 

180  PRINT  "  AS  DETERMINED  USING  GRAPHER  POLYNOMIAL  FIT  " 

190  PRINT 

200  PRINT  "  INPUT  THE  NAME  OF  THE  FILE  THAT  HAS  THE  COEFFICIENTS  FOR  THE  " 
210  PRINT  "  FOURTH  ORDER  POLYNOMIAL  FIT  WITH  ORTHOGONAL  COEFFICIENTS  " 

220  PRINT  "  IN  THE  FORMAT  AS  INDICATED  ABOVE 
230  INPUT  "",A$ 

240  OPEN  A$  FOR  INPUT  AS  #1 
250  INPUT  #1,A#,B#,C#,D#,E# 

260  INPUT  #1,A1#,A2#,A3#,A4#,A5# 

270  INPUT  #1,B1#,B2#,B3#,B4#,B5# 

280  INPUT  #1,F#,G#,H#,I#,J# 

290  INPUT  #1,A11#,A12#,A13#,A14#,A15# 

300  INPUT  #1,B11#,B12#,B13#,B14#,B15# 

310  PRINT 
320  CLS 

330  PRINT  "  THE  COEFFICIENTS  FOR  THE  FIRST  EQUATION  FOR  THE  POLYNOMIAL  FIT" 
340  PRINT  "  IN  THE  FORM  EX^4+DX^3+CX^2+BX+A  ARE  " 

350  PRINT  E#;"X^4  +  ";D#7"X^3  +  ";C#;"X^2  +  ";B#;"X  +  ";A# 

360  PRINT 

370  PRINT  "  THE  VALUES  FOR  THE  ALPHAS  IN  THE  ORTHOGONAL  FIT  FOR  THE  FIRST" 
380  PRINT  "  EQUATION  ARE  " ;A1# ;A2# ;A3# ;A4# ;A5# 

390  PRINT 

400  PRINT  "  THE  VALUES  FOR  THE  BETAS  IN  THE  ORTHOGONAL  FIT  FOR  THE  FIRST  " 
410  PRINT  "  EQUATION  ARE  " ; Bl# ; B2 # ; B3# ; B4# ; B5# 

420  PRINT 

430  PRINT  "  THE  COEFFICIENTS  FOR  THE  SECOND  EQUATION  FOR  THE  POLYNOMIAL  FIT" 
440  PRINT  "  IN  THE  FORM  EX^4+DX^3+CX^2+BX+A  ARE  " 

450  PRINT  J#;"X"'4  +  ";I#;"X^3  +  ";H#;"X^2  +  ";G#;"X  +  ";F# 

460  PRINT 

470  PRINT  "  THE  VALUES  FOR  THE  ALPHAS  IN  THE  ORTHOGONAL  FIT  FOR  THE  SECOND" 
480  PRINT  "  EQUATION  ARE  " ; All# ; A12 # ; A13 # ; A14 # ; A15# 

490  PRINT 

500  PRINT  "  THE  VALUES  FOR  THE  BETAS  IN  THE  ORTHOGONAL  FIT  FOR  THE  SECOND" 
510  PRINT  "  EQUATION  ARE  " ;B11# ;B12# ;B13# ;B14# ;B15# 

520  PRINT 

530  PRINT  "  INPUT  THE  LOWER  AND  UPPER  BOUND  FOR  THE  FREQUENCY  FOR  WHICH" 

540  PRINT  "  PROGRAM  WILL  SEARCH  FOR  A  SOLUTION 
550  INPUT  "",XL#,XR# 

555  XM#=(XL#+XR#)/2 

556  FL#=XL# 

557  FR#=XR# 

560  CLOSE  #1 
580  XL4#»E# 
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590  XL3#*D#-E#*A5#-A4#*E#-A3#*E#-A2#*E# 

600  XL2#*E#*A5#*A4#-D#*A4#+C#-B4#*E#~A3#*D#+A3#*A5#*E#+A3#*A4#*E# 

610  XL2#=XL2#-E#*B3#-A2#*D#+A2#*A5#*E#+A2#*A4#*E#+A2#*A3#*E# 

620  XL2#=XL2#-B2#*E# 

630  XL1#=A3#*A4#*D#-A3#*A4#*A5#*E#-A3#*C#+A3#*B4#*E#-B3#*D#+A5#*B3#*E# 

640  XL1#=XL1#+B#+A2#*A4#*D#-A2#*A4#*A5#*E#-A2#*C#+A2#*B4#*E#+A2#*A3#*D# 

650  XL1#=XL1#-A2#*A3#*A5#*E#-A2#*A3#*A4#*E#+A2#*B3#*E#-B2#*D# 

655  XLl #=XL1 #+A5  # * B2  #  *E #+A4  #  * B2  #  *E # 

660  XLO #=A2  # *A3  # *A4  # *A5 #  *E# -A2  #  *A3  #  *A4  #  *D#+A2  #  *A3  #  *C# 

670  XL0#=XL0#-A2#*A3#*B4#*E#+A2#*B3#*D#-A2#*A5#*B3#*E#-A2#*B# 

680  XL0#=XL0#+A#+A4#*B2#*D#-A4#*A5#*B2#*E#+B2#*B4#*E#-B2#*C# 

700  XR4#=J# 

710  XR3#=I#-J#*A15#-A14#*J#-A13#*J#-A12#*J# 

720  XR2#=J#*A15#*A14#-I#*A14#+H#-B14#*J#-A13#*I#+A13#*A15#*J#+A13#*A14#*J# 

730  XR2#=XR2#-J#*B13#-A12#*I#+A12#*A15#*J#+A12#*A14#*J#+A12#*A13#*J# 

74l  ...12#=XR2#-B12#*J# 

750  XR1#=A13#*A14#*I#-A13#*A14#*A15#*J#-A13#*H#+A13#*B14#*J#-B13#*I# 

760  XR1#=XR1#+A12#*A14#*I#-A12#*A14#*A15#*J#-A12#*H#+A12#*B14#*J#+A12#*A13#*I# 
770  XR1#«XR1#-A12#*A13#*A15#*J#-A12#*A13#*A14#*J#+A12#*B13#*J#-B12#*I# 

775  XR1#=XR1#+A15#*B12#*J#+A14#*B12#*J#+A15#*B13#*J#+G# 

780  XR0#=A12#*A13#*A14#*A15#*J#-A12#*A13#*A14#*I#+A12#*A13#*H# 

790  XR0#=XR0#-A12#*A13#*B14#*J#+A12#*B13#*I#-A12#*A15#*B13#*J#-A12#*G# 

800  XR0#=XR0#+F#+A14#*B12#*I#-A14#*A15#*B12#*J#+B12#*B14#*J#-B12#*H# 

830  PRINT 

850  PRINT  "INPUT  THE  VALUE  OF  X-MID  AND  XSCAL  FOR  THE  FIRST  " 

860  PRINT  "  EQUATION 
870  INPUT  ”",XMD1#,XSC1# 

875  PRINT 

880  PRINT  "INPUT  THE  VALUE  OF  X-MID  AND  XSCAL  FOR  THE  SECOND  " 

890  PRINT  "  EQUATION  "; 

900  INPUT  "",XMD2#,XSC2# 

905  CLS 

910  'XXLL#=(XL#-XMD1#) *XSC1# 

920  'XXRL#=(XL#-XMD2#) *XSC2# 

930  XXLM#=(XM#-XMD1#) *XSC1# 

940  XXRM#=(XM#-XMD2#) *XSC2# 

950  'XXLR#=(XR#-XMD1#) *XSC1# 

960  'XXRR#=(XR#-XMD2#) *XSC2# 

970  VAML#=  XL4#*XXLM#^4  +  XL3#*XXLM#^3  +  XL2#*XXLM#-2  +  XL1#*XXLM#  +XL0# 

980  VAMR#=  XR4#*XXRM#^4  +  XR3#*XXRM#^3  +  XR2#*XXRM#-2  +  XR1#*XXRM#  +XR0# 

990  IF  VAMR#-VAML#<0  GOTO  1020 

991  '  S#=XR4#-XL4# 

992  >  T#=XR3#-XL3# 

993  '  U#=XR2#-XL2# 

994  '  V#=XR1#-XL1# 

995  '  W#=XR0#-XL0# 

1000  XL#=  XM# 

1005  XM#=(XL#+XR#)/2 
1010  GOTO  1040 
1020  XR#=  XM# 

1030  XM#=(XL#+XR#)/2 

1040  IF  (XR#-XL#)<1E-14  GOTO  1192 

1050  IF  ABS(VAML#-VAMR#)<1E-14  GOTO  1192 

1160  GOTO  910 

1192  PRINT  "THE  SOLUTION  TO  THE  FREQUENCY  THAT  CORRESPONDS  TO  THE  INTERSECTION" 

1193  PRINT  "OF  THE  TWO  QUARTIC  EQUATIONS  IS  ";XM# 

1194  PRINT 

1195  V1#»XR4#*XXRM#^4+XR3#*XXRM#^3+XR2#*XXRM#^2+XR1#*XXRM#+XR0# 

1196  V2#=XL4#*XXLM#^4+XL3#*..XLM#^3+XL2#*XXLM#^2+XL1#*XXLM#+XL0# 

1197  PRINT  "THE  VALUE  OF  THE  NATURAL  LOG  OF  THE  AMPLITUDE  FROM  EQUATION  1  &  2" 
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1198 

1199 

1200 
1202 
1205 
1210 
1220 
1225 
1230 
1240 
1250 
1260 
1270 
1280 
1300 
1310 
1320 
1330 
1340 
1350 
1360 
1370 
1380 
1390 

1492 

1493 

1494 

1495 

1496 
1500 
1510 
1520 
1530 
1540 
1550 
1560 
1570 
1580 
1590 
1600 
1610 

1692 

1693 

1694 

1695 

1696 
1700 
1710 
1720 
1730 
1740 
1750 
1760 
1770 
1780 
1790 
1800 
1810 
1820 
1830 


PRINT  "RESPECTIVELY  ARE  ";V1#,V2# 

PRINT 

FE1#=.707*EXP(V2#) 

FE2#=  .707*EXP(V1#) 

LN2#=LOG(FE2#) 

LN#=L0G(FE1#) 

XL01#=XL0#-LN# 

XR01#=XR0#-LN2# 

XR#=XM# 

XL#=FL# 

XMl#=(XR#+XL#)/2 
XXLL#=(XL#-XMD1#) *XSC1# 

XXLM#=(XM1#-XMD1#) *XSC1# 

XXRM#=(XR#-XMD1#) *XSC1# 

VAI21#=XL4#*XXLM#^4  +  XL3#*XXLM#"3  +  XL2#*XXLM#"2  +  XL1#*XXLM#  +  XLOl# 

IF  ABS(VALM#)<1E-14  GOTO  1492 

IF  (XR#-XL#)<1E-14  GOTO  1492 

IF  ABS(VALM#)=VALM#  GOTO  1370 

XL#=XM1# 

XMl#=(XL#+XR#)/2 
GOTO  1260 
XR#=XM1# 

XMl#=(XR#+XL#)/2 
GOTO  1260 

PRINT  "THE  FREQUENCY  CORRESPONDING  TO  .707  OF  THE  PEAK  AMPLITUDE  " 
PRINT  "OF  THE  LEFT  SIDE  OF  THE  PEAK  IS  ";XM1# 

PRINT 

XL#=XM# 

XR#=FR# 

XM2#=(XM#+XR#)/2 
XXRM#=(XM2#-XMD2#) *XSC2# 

VARM#=XR4#*XXRM#^4  +  XR3#*XXRM#^3  +  XR2#*XXRM#^2  4-  XR1#*XXRM#  +  XROl# 

IF  ABS(VARM#)<1E-14  GOTO  1692 

IF  (XR#-XL#)<1E-14  GOTO  1692 

IF  ABS(VARM#)=VARM#  GOTO  1590 

XR#=XM2# 

XM2#=(XL#+XR#)/2 
GOTO  1510 
XL#=XM2# 

XM2#=(XL#+XR#)/2 
GOTO  1510 

PRINT  "THE  FREQUENCY  CORRESPONDING  TO  .707  OF  THE  PEAK  AMPLITUDE  " 
PRINT  "OF  THE  RIGHT  SIDE  OF  THE  PEAK  IS  "?XM2# 

PRINT 

LF#  =  (XM2#-XM1#)/XM# 

PRINT  "  THE  LOSS  FACTOR  FOR  THIS  DATA  IS  ";LF# 

OPEN  A$  FOR  OUTPUT  AS  #2 
PRINT  #2,A#,B#,C#,D#,E# 

PRINT  #2,A1#,A2#,A3#,A4#,A5# 

PRINT  #2,B1#,B2#,B3#,B4#,B5# 

PRINT  #2,F#,G#,H#,I#,J# 

PRINT  #2,A11#,A12#,A13#,A14#,A15# 

PRINT  #2 , Bll# , B12 # , B13 # , B14 # , B15# 

PRINT  #2,”" 

PRINT  #2, "THE  LOSS  FACTOR  FOR  THIS  DATA  IS  ";LF# 

PRINT  #2, "THE  VALUE  OF  THE  PEAK  HEIGHT  IS  ";V1# 

PRINT  #2, "AT  A  FREQUENCY  OF  ";XM# 

PRINT  #2, "FREQUENCIES  AT  .707  OF  PEAK  HEIGHT  ARE  ";XM1#;XM2# 

PRINT  #2, "THE  FILE  NAME  FOR  THIS  DATA  IS  ";A$ 

CLOSE  #2 
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1850  INPUT  "WOULD  YOU  LIKE  TO  DO  ANOTHER  CURVE  FIT  ";A$ 

1860  IF  A$="Y"  OR  A$="y"  THEN  GOTO  1880 

1870  GOTO  1900 

1880  CLS 

1890  GOTO  80 

1900  END 
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